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HMC histories for light dynamical overlap fermions tend to stay in a fixed topological sector for 
many trajectories, so that the different sectors are not sampled properly. Therefore the suitable 
summation of observables, which have been measured in separate sectors, is a major challenge. 
We explore several techniques for this issue, based on data for the chiral condensate and the 
(analogue of the) pion mass in the 2-flavour Schwinger model with dynamical overlap-hypercube 
fermions. 
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1. The 2-flavour Schwinger model 

The Schwinger model (QED2) on a Euclidean plane is characterised by the Lagrangian 

JSf(*,V ) A M ) = *( J r)[r M (id M + ^ M )+m]V(jc) + iF MV (jc)F M v(jc) . (1.1) 

Analytic results are obtained in the bosonised form, e.g. for the chiral condensate £ = — ( , J n F) and 
the "meson" masses (iso-triplet and iso-singlet), with two degenerate flavours of mass m <C g, 



E(m) =0.388... (mg 2 ) 1 ^ 3 § , M % = 2.008 . . . (ra 2 g) 1//3 § , M, = y^+M2 g . (1.2) 

Here we consider a lattice formulation with compact link variables U^^ G t/(l) and the pla- 
quette gauge action. For the fermions we use the overlap-hypercube Dirac operator |]] 



A,viiF - (1 -| )D ( S BP + m, D£^ F = l + (A 1 F-l)/V( fl HF- 1 "Aif- 1 > • (1-3) 



Dhf is a truncated perfect hypercube fermion operator, which is y 5 -Hermitian, D HF = /sDhf75 ■ 
Due to the use of the overlap formula in eq. ( |L~3| ), D^ HF solves the Ginsparg-Wilson relation in 



its simplest form [Q] , {D^jp, 75} = 2D^ HF y5D^ HF , which implies a lattice modified, exact chiral 
symmetry It reproduces the axial anomaly correctly in all topological sectors [p. 

In the free case Dhf is discussed in Refs. |[] (we refer to the version denoted as CO-HF). It 
is approximately chiral already, hence it changes only little in the transition to the overlap operator, 

£>hf~£> ovH f, (L4) 

in contrast to the Wilson kernel, which is used in the standard overlap operator Further virtues 
of Dhf, which are based on the renormalisation group construction, like improved scaling and 
approximate rotation symmetry are essentially inherited by D ov hf- Also long-range couplings are 



only turned on slightly, again due to the similarity (|1.4|), which strongly improves the degree of 
locality compared to the standard overlap operator. As a related property, the condition number of 
the operator in the inverse square root of eq. ( |1.3[ ) is strongly reduced. 

The form of D ov hf[^] interacting through U(l) gauge fields was introduced in Ref. [Q]. All 
the above virtues were tested and confirmed extensively for the 2-flavour Schwinger model. In that 
case, the configurations were generated quenched, but their contributions to the measurements were 
re-weighted with the fermion determinant [Q, Also in quenched QCD, a drastically improved 
locality and approximate rotation symmetry have been confirmed for D ov hf[^] iJTQj] . 

2. Simulation with dynamical overlap fermions 

Hybrid Monte Carlo (HMC) simulations with chiral fermions are tedious and still at an early 



stage. However, in addition to the virtues listed above, the property (1.4) also facilitates HMC 
simulations: a low polynomial in Dhf can be used for the HMC force and for short trajectories 
a useful acceptance rate persists [ p| , |12| ] (in the extreme case of using directly Dhf, however, its 



volume dependence was considered unsatisfactory [13]) 
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We performed HMC simulations with dynamical overlap-hypercube fermions at /3 = l/g 2 = 
ontxL lattices with two fermion flavours of mass m, where L and m take the following values: 

L= 16 : m = 0.01 , 0.03 , 0.06 , 0.09 , 0.12 , 0.18 , 0.24 . 



m 



0.01 : L = 16 , 20 , 24 , 28 , 32 



We did not study a continuum extrapolation, but the lattices are fine (plaquette values « 0.9), so 
we rely on small lattice artifacts (moreover 0(a) scaling artifacts are ruled out since we are using 
Ginsparg-Wilson fermions). 

Our algorithm fulfils the conceptual conditions like detailed balance and area conservation 



[|11[]. The essential practical question is the acceptance rate in the Metropolis step at the end of the 
HMC trajectories (it employs D ov hf to precision 10~ 16 ). For L = 16 with trajectory length £ = 1/8 
we obtained acceptance rates in the range 0.3 .. .0.5, with only a mild dependence on the fermion 
mass [11]. In Figure [T] we show the dependence on the lattice size L at very light mass m = 0.01. 
The trajectory length is decreased for L= 16, 20, 24, 28, 32to^ = 0.125, 0.0625, 0.05, 0.04, 0.03, 
respectively. This keeps the acceptance rates in the same range, and large statistics allows for a se- 
lection of well decorrelated configurations. 

A major challenge for the simulation 
with light dynamical overlap fermions — 
in addition to the huge computation time 
request in QCD — is that the histories at 
weak gauge coupling perform only very few 
topological transitions (which require a tem- 
porary deformation to rough configurations). 
We refer to the clean definition, which iden- 
tifies the topological charge with the fermionic 
index v [^J. In all our considerations only 
the absolute value |v| matters. 



S 

Q. 

CD 



0.6 



0.5 



0.4 



0.3 - 



0.2 



~i 1 1 

m = 0.01 



15 



20 



25 
L 



30 



35 



3. The chiral condensate E 

For a configuration in volume V 
the chiral condensate is given by 



Figure 1: The acceptance rate of the HMC trajectories at 
m = 0.01 as a function of the lattice size L. Note that the 
trajectory length is reduced for increasing volume. 
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(3.1) 



A,- are the Dirac eigenvalues mapped stereographically onto the imaginary axis, A,- — ► A,/(l — A,/2), 
and the sum runs over all of them. The last term defines the quantity £i v i . 

In Ref. [[14]] we discussed the determination of £ based on the lowest non-zero eigenvalues and 
Random Matrix Theory (RMT). Since E(m = 0) = 0, this is not the setting that RMT usually refers 
to, so we are probing terra incognita. In fact we observed microscopic eigenvalue densities, which 
are not described by any known RMT formula. Nevertheless, the RMT interpretation of the ratio 
(Ai)v=o/(Ai)| v |=i yields results for E, which agree well with the predicted value (T2) over abroad 
parameter range. ((. . . )i v i denotes an expectation value restricted to configurations of charge |v|.) 
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Here we address the direct measurement of £ based on the full Dirac spectrum. Our results 
are shown in Figure At small masses m, the zero mode contribution strongly dominates for 
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Figure 2: The contributions of the non-zero modes to the chiral condensate, £| V |, measured in the topological 
sectors ±v. We show the mass dependence at L = 16 (left) and the lattice size dependence at m = 0.01 (right). 

topologically charged configurations, so that the term ei v i is a minor correction. As a generic 
property of stochastic Hermitian operators (such as 75D ov hf)> zero modes repel low lying non-zero 
modes, which implies the hierarchy 

£ > £\ > s 2 . . • (3.2) 

at fixed m and L. This is observed consistently in Figure || (on the left). As we increase the volume 
at fixed m, more eigenvalues are accumulated near zero, and we infer 

e*(Vi) > ei(y 2 ) for Vx>V 2 . (3.3) 
Also that property is observed consistently from our data, as Figure ^| (on the right) shows. 

3.1 Method 1 : Gaussian summation 

Our first method to sum over the topologies relies on the assumption that the probability dis- 
tribution of the charges v is Gaussian. Indeed general experience shows that possible deviations 
from a Gaussian distribution tend to be small. Here our assumption implies 



I = £ P(|V|)L| V | , where p(|v| 



exp{-v 2 /(2V Z? )} 
L v exp{-v 2 /(2^)} ' 



-(yy)|v| anc ^ Xt = (v 2 ) /V is the topological susceptibility. 



(3.4) 



Assume we have measured results for Lq...Lq. This constrains the values at all higher charges 

de- 



based on eq. (3.2), 



]V|_ < £ |V| 

mV ' v ' mV 



(3.5) 



For L = 24 and 28 we also have gaps at low | v | , where we insert the obvious bounds based on eq. 



(3.3). Thus we obtain — with remarkably mild uncertainties — the full set of Ei v i to be inserted 



in eq. (3.4). The only unknown parameter is the susceptibility % t , which we tune so that the sum 



4 



Topological Summation with Dynamical Overlap Fermions 



Wolfgang Bietenholz 



0.002 



0.0015 



0.001 



0.0005 



I I / I ,. I 
L- 16 / / 










I " 






0.0067 m 1 2 




/ conjecture (1/ Xt Q = 0) -- 
(X-.u-Wl) 







0.02 



0.04 0.06 
m 



0.08 0.1 



0.0007 
0.0006 
0.0005 
0.0004 
0.0003 
0.0002 
0.0001 




I 


i 


E (1) m/N f 

i 


m = 0.01 - 

i 



400 



800 



1200 



V 



Figure 3: Results for the topological susceptibility % t for light fermions, obtained by assuming a Gaussian 
distribution of the topological charges and matching the theoretical value of E. The plot on the left suggests a 
behaviour % t °< %fm at small masses. The plot on the right shows a nice convergence to the value conjectured 
in Ref. Jl6|] (in infinite volume), from where we also adapted the quenched susceptibility % t ,q- 



reproduces the theoretical Z value (1.2). Thus we obtain the results for % t shown in Figure |3| For 
very light fermions they suggest %, oc ^/in (in a fixed volume), see plot on the left. 

Alternative results (with quenched configurations and re- weighting) were given in Ref. [|15|] . 
On the theoretical side, Ref. [jlfj conjectured for Nf degenerate flavours in the large volume limit 



1 

Xt 



N f 1 
+ — 
m Xu 



:=L N i(m = 0)~0.16g , 



(3.6) 



where % t q = Xt( m —> °°) 1S the quenched value. As we increase the volume at m = 0.01, our results 
converge to the vicinity of this prediction, see Figure |3| on the right. 



3.2 Method 2 : Summation formula 

The first method used the theoret- 
ically predicted Z value as an input to 
extract % t . Now we try to evaluate £ it- 
self from our data. We apply an approx- 
imation formula [O], in exact analogy 
to a formula for the pion mass derived 



in Ref. [17], 



0.09 



0.08 



0.07 



0.06 



0.05 



£ -77 + V 2 77T, 



V 



V 2 



Xt 



8=4 
Xt 

(3.7) 



0.04 



n 1 1 1 r 

T theory (V infinite) 



m = 0.01 



j_ 



200 400 



600 
V 



800 1000 1200 



The derivation involves a Fourier trans- Figure 4: The chiral condensate E at m = 0.01 as a function 
form between V and the vacuum angle of the volume, obtained from measurements of various Ei v j 



6, which is treated by the stationary phase values and the approximate summation formula (f3/7|). 
approximation. 

Formula (3.7) involves three unknown parameters. Two of them, £ and Xt = A/B are of physical 
interest (unlike a). We capture them if we manage to determine A and B. 
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Figure 5: M n ' and as functions of the fermion mass m at L = 16, measured in the sectors |v| = and 
1. For comparison we show the theoretical prediction (summed over all topologies). 

• At fixed m and V we can extract B, e.g. from Eo and Ei. 

• At fixed m in two volumes V\ ^ V%, we can evaluate A, e.g. from Eo(Vi) and Eq^)- 
Ei v i values in three volumes determine all three parameters. 



For a combined approach along these lines [12], we arrived at the results for E shown in Figure 
|J As in Figure || (on the right) we see for increasing volume a flow towards the value (1.2), which 
was predicted theoretically (in infinite volume). 

3.3 Method 3 : Density correlation 



A third method has been suggested in Refs. [18]. The derivation is similar to Ref. [[170, 1 an d 
it leads to a formula for finite size effects in the correlation of the topological charge density p , 



<P(x)p(y)>[ 



V 



Xt 



C4 

2XtV 



+ 0(V~ 



(3.8) 



(The kurtosis C4 = (3(v 2 ) 2 — (v ))/V is a measure for the deviation from a Gaussian charge dis- 
tribution.) The formula refers to the asymptotic behaviour at large |jc — y\. This requires a large 
volume, and unfortunately this also means that the signal to extract % t is likely to be too small for a 
conclusive measurement. Later Ref. [|l^] considered instead the r\ '-correlator of the pseudo-scalar 
density, which obeys the analogous formula, and where a sensible signal was found, even in QCD 



with dynamical overlap quark. Also in our project this method is under investigation [12]. 



4. Meson masses 

We measured the "meson masses" based on current correlators, which is most efficient in this 
model [Q|, and we show the results for L = 16 in Figure ||. At least for m = 0.01 we are clearly in 
the £ -regime. For increasing m the topological distinction shrinks. The masses measured in | v | = 
and 1 are compared to the predictions of Refs. [jl[ ||] (which refer to m <C g — 0.45). 



At last we apply Method 2 to the pion mass, as it was originally intended [|17p. We fix m = 0.01, 
= 1 and employ the formula (3.7), M n \ w M % — A/V + B/V 2 . We have results from three 

'in both cases the approximation should be best for small |v| , so that (v 2 ) 2> |v|. 
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volumes, which we insert and solve for M % , 
{ M^ =16) = 0.276(4) , M { n L = 20) = 0.214(4) , m[ L = 32) = 0.135(4) } =^ M n = 0.078(8) . (4.1) 



This is indeed compatible with the theoretical value in eq. (1.2), M n = 0.0713 . . . Considering the 



large masses that we started from, this agreement is an impressive success of this method. 
5. Conclusions 

Based on our experience, it appears possible — at least in some cases — to derive observables, 
which are properly summed over all topologies, even if only measurements in a few fixed sectors 



are available. This property (which agrees with Ref. [19]) is crucial for the future of dynamical 
overlap fermion simulations, in particular in the £ -regime, where topology is essential pfj. The 
exact results and their reliability depends on subtleties and needs to be explored further — the 
Schwinger model is ideal for such tests before large-scale QCD applications. 
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going project, and S. Diirr for helpful comments. The simulations were performed on the clusters 
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